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Studies have shown that emission of gravitational waves drives an instability in the r-modes 
of young rapidly rotating neutron stars carrying away most of the angular momentum through 
gravitational wave emission in the first year or so after their formation. Magnetic field plays a 
crucial role in the evolution of these r-modes and hence the evolution of the neutron star itself. An 
attempt is made here to investigate the role of magnetic field in the evolution of r-mode instability 
and detectibility of gravitational waves emitted by a newly born, hot and rapidly and differentially 
rotating neutron star. It is found that magnetic fields tend to suppress the r-mode amplitude. The 
signal-to-noise ratio analysis shows that gravitational waves emitted from the r-mode instability 
from neutron stars with magnetic fields upto the order of 10 14 gauss may be detectable by the 
Advanced LIGO at 20 Mpc. 

PACS numbers: 04.30.Tv, 95.85.Sz 



I. INTRODUCTION 

Recent developments in the kilometer— scale laser in- 
terfcrometric gravitational— wave detectors LIGO and 
Virgo, or their advanced versions, operating in the fre- 
quency range of 10 and 10 Hz usher a new era in the 
gravitational-wave astronomy, enhancing the prospects 
of detecting gravitational waves from a variety of astro- 
physical sources. 

The r-modes are primarily non-radial velocity pertur- 
bations in stars. These perturbations are of the order of 
the angular velocity f2 of stars' rotation and entail density 
perturbations of order Q 2 . The characteristic frequency 
of oscillation is also comparable to fl. 

The interesting fact about the r-modes is that they 
are driven unstable by gravitational radiation reaction 
in all perfect-fluid rotating stars for a large range of rel- 
evant angular velocities. This property was discovered 
by AnderssonQ and confirmed theoretically by Fried- 
man and MorsinkQ. So, unlike other modes, r-modes 
are always retrograde in the star's corotating frame and 
progradc in the inertial frame, i.e., the sign of the r-mode 
frequencies is always opposite in the two frames. In other 
words the Chandrasekhar-Friedman-Schutz (CFS) insta- 
bility always occurs for the r-modes. 

The effects of bulk and shear viscosity in neutron star 
matters are to dampen the driving effects of gravitational 
radiation. It is found that in case of neutron star matters, 
the effect of gravitational radiation driven instability is 
stronger than the damping effects of viscosity P, Q ■ 

A phenomenological model for evolution of the r-mode 
instability was given by Owen et. aZ.[j| to study the de- 
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tcctibility of gravitational waves emitted by a newly born, 
hot, rapidly-rotating neutron star. They found that such 
waves could be detected by the enhanced version of LIGO 
if sources were located at distances up to 20 Mpc from 
Earth. However, a deeper understanding of this issue re- 
quires taking into account of nonlinear effects in the evo- 
lution of r-mode instability. Arras et al. Q studied the 
nonlinear mode interaction to conclude that enhanced 
laser intcrfcromctric detectors could detect gravitational 
radiation from r-modes provided the sources are located 
within a distance of 200kpc. That r-modes can induce 
a drift of fluid elements along azimuthal directions was 
suggested by Rezzolla et al. who derived an approxi- 
mate analytical expression for these drifts and these were 
numerically confirmed in both gen eral relativistic Q and 
Newtonian Hydrodynamics [Tol 1 1 1| regimes. 

Recently, Sa [12| has found an exact r-mode solu- 
tion within a nonlinear theory up to second order in the 
mode's amplitude which describes differential rotation of 
pure kinematic nature producing large scale drifts along 
stellar latitudes. This solution is then used to investigate 
the influence of differential rotation in the evolution of r- 
mode instability of a newly born, hot, rapidly rotating 
neutron star [13[ and the detectibility of gravitational 
waves emitted by such a system by Advanced laser in- 
terometric detector LIGO [l4j . 

A newly born hot neutron star is beset with high spin 
and coherent, large-scale hot plasma fluid currents with 
high electrical conductivity. In such neutron star mat- 
ter, large scale magnetic fields can be produced. It will 
thus be interesting to investigate the influence of mag- 
netic field in the evolution of newly born hot neutron 
star, in particular its influence on the spin evolution 
through r-mode instability. Spruit (l5j first investigated 
the generation of a toroidal magnetic field as a result of 
gravitational wave emission induced differential motion 
in the star. The huge magnetic field produced along with 



the buoyancy instability associated with them have been 
used for the production of a gamma-ray burst model. Ho 
and Lai [l(| introduced magnetic braking originated from 
the magnetic dipole radiation in the spin evolution of a 
unstable neutron star and estimated the growth rate of 
the instability as a result of fast magnetosonic and Alfven 
wave emission in the star's magnetosphere. 

The present work attempts to make phcnomenological 
descirption of the nonlinear evolution of the r-modes in a 
newly born hot neutron star up to second order in mode's 
amplitude by extending the work of Sa and investigate 
the influence of the magnetic field on the mode amplitude 
and the kinematic differential rotation. 

In section II we describe briefly the general properties 
of the r-mode, section III details evolution of r-modes 
in a canonical neutron star with or without invoking the 
differential rotation and magnetic fields. In section IV 
we calculate the waveform of gravitational waves emitted 
during r-modc evolution. The results are analysed in 
section V. We summarise our results in section VI. 



II. THE R-MODE INSTABILITY: GENERAL 
CONSIDERATION 



the energy E in the rotating frame and the time deriva- 
tive of E. The energy E can be expressed as a quadratic 
functional of fluid velocity perturbations as: 
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We can decompose 1 /r as 
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where 1/tgr,1/ts and 1/tb are the contributions due to 
gravitational radiation emission, shear viscosity and bulk 
viscosity, respectively. Expressions for these individual 
contributions for the r-modes are given by [l[ 
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The Eulerian velocity perturbation in the rotating 
Newtonian stars in case of r-modes are given by 
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Sv = RQf(r/R) Y& 
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where R and f2 are the radius and angular velocity of the 
unperturbed star, f(r/R) is an arbitrary dimensionlcss 
function, and is the magnetic type vector spherical 
harmonic defined by 
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For barotropic models, the radial dependence f(r/R) 
equals a(r/R) 1 , where a is an arbitrary constant [18l ]. 
These modes exist with velocity perturbations as given 
by Eq. (1) if I = m, and the mode frequencies are 
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Since the Coriolis force serves the restoring force for these 
oscillations, the frequencies of these modes are low com- 
pared to the / and p -modes in slowly rotating stars. 
These expressions for Sv and u) are the lowest order in an 
expansion in terms of the angular velocity Q. The exact 
expressions contain additional terms of the order Q 3 . 

The form of the r-mode solution has the time depen- 
dence e ™ t- */ T as a consequence of ordinary hydrodynam- 
ics and the influence of the various dissipative processes. 
The real part is the frequency of these modes, (J, given 
in (3), while the imainary part 1/r is determined by the 
effects of gravitational radiation, viscosity, etc. It is pos- 
sible to evaluate 1 /r by computing the time derivative of 
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where dp is proportional to O 2 and hence is small (i.e., 
higher order in J7) compared to Sv in slowly rotating 
stars. 

The expression for 1/rg in Equation (J9j> is only ap- 
proximate. The exact expression should contain the La- 
grangian density perturbation Ap in place of the Eulerian 
perturbation Sp. The bulk viscosity [Eq. pTjl ] is a very 
strong function of the temperature, being proportional 
to T 6 . 

Owenet. al.(1998) @ have evaluated these expressions 
for the imaginary parts of the frequency for a "typical" 
neutron star model with a polytropic equation of state : 
p = kp 2 , with k chosen so that a 1.4M© model has the 
radius 12.53 km by using the ex pre ssions for the viscosity 
of the hot neutron star matter UM, uOl : 
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where all quantities are expressed in cgs units. 

It will be useful to define a timescale associated with 
the viscous dissipation 1 fry = 1/rg + 1/tb- The viscous 
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timescale ry and the gravitational timescale tqr can be 
expressed in terms of the temperature and angular veloc- 
ity dependences as 
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where ts = 2.52 x 10 8 seconds, tb = 6.99 x 10 8 seconds 
and tqr = —3.26 seconds for I = 2 mode 



III. EVOLUTION OF THE R-MODES IN 
NEUTRON STARS 

A. In the absence of differential rotation and 
magnetic field 

In the consideration of the evolution of the r-modes 
in neutron star environments, it is assumed that initially 
the mode will be a small perturbation that is described 
adequately by the linear analysis. However, as the mode 
grows, non-linear hydrodynamic effects become impor- 
tant and eventualy dominate the dynamics. In the ab- 
sence of the mathematical tool to exactly follow this non- 
linear phase of the evolution, Owen et. al. Q first pro- 
posed a simple linearised treatment that provides the ba- 
sic features of the exact evolution. Essentially the neu- 
tron star is treated as a simple system having only two 
degrees of freedom: the uniformly rotating equilibrium 
state parametrized by its angular velocity f2, and the r- 
mode is parametrized by its amplitude a. 

The canonical angular momentum, i.e., the angular 
momentum associated with the r-mode, can be expressed 
in terms of the velocity perturbation Sv by 
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For the I = 2 r-mode, this expression at the lowest order 
in ft reduces to 
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where J is given by 
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For polytropic models the dimensionless constant J = 
1.635 x 10 -2 . The moment of inertia / can be written as 



/ = I MR 2 



where / is given by 
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and for r = 2 polytropic models, / ~ 0.261. 

The angular momentum associated with the mean rigid 
rotation primarily increases through a transfer from the 
mode due to viscosity through the emission of gravita- 
tional radiation. Thus, 
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The rate of change of canonical momentum of the mode 
can increase through gravitational radiation and decrease 
by transferring angular momentum to the star through 
viscosity. Consideration of this fact leads to another evo- 
lution equation for the r-modes. This was realised by 
Ho and Lai [l6[ that it is the evolution of the canonical 
momentum and not the energy evolution as in Owen et. 
al.Q, that is responsible for the correct description of the 
evolution of r-mode. Thus, 
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Substituting equations ([T5"]) through (|T5|) into equation 
(PI) yields 
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B. r-mode evolution of differentially rotating 
neutron stars but no magnetic field 

Sa et. al. [l2|, [l3[ while attempting to find out the 
nonlinear correction (in the second order in the mode's 
amplitude a), to the equations of Owen et. al, found that 
the leading non-linear effect can be interpreted as dif- 
ferential rotation only if one includes a correction term 
to the expression of the canonical angular momentum. 
After a detailed nonlinear calculation, they found that 
the corrected version of the canonical angular momen- 
tum term (called as physical angular momentum of the 
r-mode) is given by [l3| 
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where Q, R and Al denote the angular velocity, the radius 
and the mass of the star, respectively so that the total 
angular momentum of r-mode is 



(25) 



and K is a parameter accounting for the differential ro- 
tation in the r-mode evolution. 

Incorporation of the above correction leads to the r- 
mode evolution equation for a differentially rotating star 
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It is to be noted that the evolution equations for a rigidly 
rotating star (similar to Q and also [16[ with the suppres- 
sion of the magnetic field term) can be recovered from the 
above equations by setting K = —2. Another advantage 
with this set of evolution equation is that they account 
for the entire history of r-mode evolution - one need not 
monitor the saturation value of a and then change the 
governing equation as had to be done in the case of either 
Owen et. al or of Ho and Lai. 



C. r-mode evolution of differentially rotating 
neutron star in the presence of magnetic field 

As the neutron star is beset with a very strong 
magnetic field, of the order of 10 14-16 Gauss , the role 
of such a huge magnetic field cannot be ruled out in 
governing the evolution of r-modc. This is quite possible 
that evaluation of different other quantity like the 
amplitude of gravitational wave that is emitted by such 
a neutron star during r-mode evolution period and the 
amount of spin down and spin down period etc. will be 
overestimated if the magnetic field is not accounted for 
during ?*-mode evolution of a neutron star. Ho and Lai 
incorporated the effect of magnetic field in a very simple 
way by considering the magnetic braking effect during 
the spin down of the star. Though a detail analysis of 
the role of magnetic field is itself a formidable task, in 
this paper we incorporate the magnetic field into the 
evolution equation of Sa et. al. in the same line as done 
by Ho and Lai and analyse its effect in the differential 
rotation of the star during r-mode evolution. We further 
calculate the amount of gravitational radiation to be 
emitted by a neutron star in the newtonian regime with 
a polytropic equation of state and for different strengths 
of magnetic fields and differential rotation parameter 
to have a qualitative estimate of the role played by the 
magnetic field and differential rotation simultaneously 
in the gravitational wave emission. We further look into 
the gravitational wave detectibility aspect of such a star 



in the eyes of LIGO I, LIGO II and VIRGO. 

In the I = 2 r-mode instability of second order non- 
linear theory, two two crucial parameters are the angular 
velocity fl and the mode amplitude a. The physical an- 
gular momentum of the second order perturbation, which 
includes the differential rotation in a typical neutron star 
is given by 



j( 2 ) J = ia 2 fi(4A + 5)JAIR 2 



(28) 



This expression assumes the star's mass density p and 
pressure p to be related by a polytropic equation of state 
p = kp 2 and the constant k is chosen such that M = 
1.4Mq and R = 12.53 kilometres. The constant K in 
equation (j2"5)l is to be fixed by initial data, giving the 
initial amount of differential rotation associated with the 
r-mode (l2j . 

The total angular momentum is the sum of the bulk 
angular momentum and the mode angular momentum, 
or the angular momentum of oscillation, is given by 
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where I = (8tt/3) /„ pr 2 dr = I AIR 2 , {I = 0.261) is 
the moment of inertia of the unperturbed star. 

The angular momentum of oscillation S^- 2 \l increases 
mainly due to the emission of gravitational wave and the 
dissipative effects of bulk viscosity and shear viscosity. 
Taking into consideration of these facts we find 
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On the other hand, the total angular momentum 
J(a, Q) decreases due to emission of gravitational wave 
and Magnetic braking in the case of a star beset with 
magnetic field of strength B. This gives 
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The magnetic braking term arises from the radiation 
of rotational energy from the rotating magnetic neutron 
star over a time scale tm due to the intervention of mag- 
netic field and leads to an equivalent magnetic braking 
torque of magnitude —10,/tm- For a simple magnetic 
dipole model of neutron star, we have plj 
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Thus the consideration of decrease of total angular mo- 
mentum due to emission of gravitational wave and mag- 
netic braking, leads to the equation 
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Using equation (|29[) . equation (|33|) reduces to 
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where 



On simplification, 
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On the other hand, consideration of the increase of 
5^ J(a, f2) due to gravitational radiation instability and 
decrease due to bulk and shear viscosity, leads to 
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which, when used in equation ()34[) . gives 
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Equations (|36|) and (|37j) arc the governing equations 
for evolution of r-mode instability in the presence of both 
differential rotation and magnetic field. K is a constant 
which is to be fixed by the choice of initial data, giving 
the initial amount of differential rotation associated with 
the r-mode. The value of K is chosen to lie in the interval 
—5/4 < K <gi fO 17 . The upper limit for K results from 
the imposition of the condition that the initial value of 
the angular momentum of the r-mode is much smaller 
than the angular momentum of the unperturbed star, 
i.e., 5& < Ift , which gives K <C 10 17 , if we choose ao = 
f0~ 8 . The lower limit for K results from the fact that 
we want to avoid the total angular momentum of the star 
becoming a negative quantity. This will further ensure 
that amplitude saturation is governed by the equations 
themselves unlike the case in [5j where the saturation of 
the r-mode amplitude was put by hand. 

The equations (f3l)]) and (|3"T)) are numerically solved for 
a range of values of initial differential rotation and mag- 
netic field strengths. Convergence of the numerical so- 
lutions are checked by changing the values of the step- 
size of iterations and found to be numerically convergent. 
For the purpose of illustrating the behaviour of nonlinear 



growth, we have chosen the values of differential rotation 
parameter K and the magnetic field strengths as the fol- 
lowing: 



K 



Bi 



10 100 1000 10 4 10 5 10 6 10 7 10 s 10 M 10 iu 10 11 10 



,9 mlO mil inl2 



2 4 6 8 10 15 20 40 60 80 100 



TABLE I: Table of different values of initial differential rota- 
tion parameter K and the magnetic field B\± (xlO 14 Gauss) 
chosen for the numerical solution of Equations (|36|) and (|37|) . 

B14 is the value of the magnetic field strength in units 
of 10 14 Gauss. For different combinations of the values 
of K and £?i4 chosen above, we have numerically solved 
the equations and analysed the results for the growth of 
mode amplitude, evolution of angular velocity, gravita- 
tional wave amplitude and the detectibility of emitted 
gravitational wave signals. 

It is further necessary to specify how the temperature 
of the star evolves with time in order to make the model 
of the evolutionnof the r- modes complete. We adopt here 
the standard description of the cooling of hot young neu- 
tron stars due to emission of neutrino via the modified 
URCA process. The temperature during this phase falls 
quickly by a simple power law cooling formula [2lj|. 
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where Tj is the initial temperature of the neutron star, 
and t c is a parameter that characterises the cooling rate. 
For the modified URCA process, r~ 1 year. A typical 
value for the initail temperature is Tj w 10 11 K. This 
equation can now be inserted into the evolution equations 
for a and Q, to provide explicit differential equations fo 
rthe time evolution of the angular velocity of the star and 
the amplitude of the mode. 
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FIG. 1: Growth of r-mode amplitude for different values of 
the differential rotation parameter K, and magnetic field B=0. 



IV. GRAVITATIONAL WAVEFORM FROM 
R-MODE INSTABILITY IN THE PRESENCE OF 
DIFFERENTIAL ROTATION AND MAGNETIC 
FIELD AND THEIR DETECTIBILITY 

The frequency-domain gravitational waveform 



Hf) = 



a 2nift 



h(t)dt 



(39) 



is determined completely by the assumption that the an- 
gular meomentum radiated as gravitational waves comes 
directly from (the angular momentum of) the star. This 
assumption is expected to be satisfied during the non- 
linear saturated phase of the evolution, but not during 
the early evolution when the mode is growing exponen- 
tially. This follows that the rate of the frequency of the 
radiation evolves with time in the later phase of the evo- 
lution. 

In the stationary phase approximation, the gravita- 
tional wave strain h(t) is related to its Fourier transform 
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FIG. 2: Growth of r-mode amplitude for different values of the 
differential rotation parameter K, and magnetic field B14 = 
10. 
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FIG. 3: Growth of r-mode amplitude with time for differ- 
ent values of K, the differential rotation parameter and with 
magnetic field B14 = 20. 



Hf) by 



W)\ 2 = \h(f)\ s 



(40) 



For the I = 2 of the r-mode which we discuss here, the 
mode frequency is w = |fi, or / = jj^fl, where / is the 
frequency of the emitted gravitational waves measured in 
Hz. 

The rate of angular momentum loss due to the r-mode 
can be coupled to the emission of gravitational waves in 
the standard current multipole expansion for I = m = 2 
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FIG. 4: Growth of r-mode amplitude with time for differ- 
ent values of K, the differential rotation parameter and with 
magnetic field B14 — 100. 
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FIG. 5: Growth of r-mode amplitude with time for different 
values of magnetic field and with differential rotation param- 
eter K=0. 
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FIG. 6: Growth of r-mode amplitude for different B and K 
10 6 . 
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FIG. 7: Growth of r-mode amplitude for different B and K 
10 12 . 



Again, the angular momentum loss couples to the emis- 
sion of gravitational waves as[22| . 
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where h + (t) and h x (t) are the strain amplitudes for the 
two polarisations of the gravitational wave, D is the dis- 
tance to the source and < • ■ ■ > denotes an average over 
possible orientations of the source with respect to the ob- 
server. Combining the Equations (1421) and (|43|) we arrive 
at 
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In the stationary phase approximation (which is al- 
ways valid for a secular instability), this expression for 
the gravitational wave is written in terms of the Fourier 
transformed quantity as 



<h 2 + (f) + h 2 x (f) >-- 
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can be calculated from the expres- 
sion for ^ given by the Equation (|37p and using Equa- 



The quantity 



tions (jnjl, UM and (E2J) we get 



Here fl t = fi/f^o with D,q = V^pG, p being the av- 
erage density of a typical neutron star. Also, since the 
viscous timescale tv compares very long (of the order of 
few years) with respect to either tqr or tm (of the order 
of few hours) , we have neglected the viscosity term from 
the above equations. 
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With the help of Equation (|46|) . it is now easy to rewrite 
Equation (|45|) as 
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The measured value of 



Kf) 



depends on the loca- 



tion of the source in the observer's sky and the polarisa- 
tion angle with respect to the interferometer. Since the 
gravitational wave interferometers arc omnidirectional, 
we must have to effect an averaging over these angles to 
have a correct description of the measured strain values. 
Thus, averaged over these angles, the measured strain in 
the Fourier domain is given by, 



Also, the average over the source location in three- 
dimensional space should also be taken into account 
because of the fact that spatial average weights more 
strongly those orientations giving out stronger signals 
and the effect of this averaging enhances < \h(f)\ 2 > 
further by about 1.5 (23| . Combining these results, the 
average value of h(f) generated in a neutron star like 
source reduces to 3/10-th of that given in Equation (|4T|) 
and we have, 



<K/)I 2 >=5<IM/)I 2 +<IM/)| 2 > (48) 



~ 64 J^GM&Ja^ r 
u ' 75 V 3 c 5 D * 



8 (K + 2)a 2 Qn* 

3 TGR 



B 14 
TM 



-1/2 



(49) 



r 



For a typical neutron star of mass M = lTM©, and 
radius R = 12.53 km emitting an average value of h at a 
distance of D — 20 Mpc is then given by, 



function of time can be found using Equation ([4"4")) . When 
averaged over source and detector orientation, it is given 
by 



h(f) = 3.3 x 10~ 26 



n t (K + 2) 



B 2 
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6.74 x 10 4 a 2 f2? 



-1/2 



(50) 



The characteristic gravitational wave amplitude as a 



|„ (1 )|=6.6xl0-"a(<)(^) (^), (51) 

where D is the distance to the source in 20 Mpc and 
ilo = y/irpG. We have considered the evolution for a 
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period of one year from the inception of the instability. 
Figures (fT5|) through (f2~Tj) depict the evolution of abso- 
lute value of gravitational wave strain during r-mode in- 
stability for various values of magnetic fields B and the 
differential rotation parameter K. During the early part 
of the evolution, that is, upto about 12 minutes in case 
of K=0 and about 4 minutes for K=10 12 , \h(t) \ grows ex- 
ponentially while during the later part, extending upto 1 
year, it decreases slowly. At its maximum, h(t) assumes 
a value of 5 x 1CT 25 for K=0 and 1 x 1CT 30 for K=10 12 , 
both for Bi4=0. The cases with other values of K and 
-B14 also behave similarly. 

It is worth mentioning here that in Ref. [f| an expres- 
sion for the frequency-domain gravitational wave ampli- 
tude was derived based on the assumption that 
dJ/df oc 7, where J is the total angular momentum of 
the star, / is the frequency of the emitted gravitational 
wave and I is the moment of inertia of the unperturbed 
star. Since within the model of Ref. Q, the condition 
dJ/df oc I only applies during the second stage of evolu- 
tion, the expression obtained there for \h(f)\ is valid just 
for this stage of evolution. Within our consideration it 
can be shown that the condition dJ/df oc I easily holds, 
not only during the seend stage of evolution, but also 
during the first stage. 

V. ANALYSIS OF RESULTS 

First, we choose to consider the duration of evolution 
of I = 2 r-modc instability to be a period of 1 year from 
the inception of the instability. During this period, the 
dissipative effects of bulk viscosity tb is hardly dominant 
to damp the mode. Temperature during this period falls 
from about 10 11 K to 10 9 K (Equation (|38jl). We wish 
to stop the evolution beyond 1 year since temperature 
achieved beyond 1 year of evolution, supcrfluid effects 
are expected to become important, rende ring invalid our 
assumptions about the viscous timcscale [24j. 

The amplitude of the r-mode at the inception of the 
instability is taken to be as small as 10 -8 . The problem of 
solving the Equations (|36|) and (|37| with no contribution 
from tb then is amenable to straight numerical solution 
by Runge-Kutta 4-th order method. 

The curves for growth of r-mode amplitude a, the evo- 
lution of angular velocity VL and the gravitational wave 
strain amplitude h, suggest that the entire r-mode in- 
stability growth proceeds in two stages - the first stage 
and the second stage. During the first stage, the early 
stage, which lasts about few hundred seconds, the mode 
amplitude and gravitational wave strain amplitude expo- 
nentially rise and the angular velocity maintains more or 
less a constant value. During the second stage, the later 
stage after the first stage, lasting for the rest period of 1 
year, the angular velocity falls abruptly. The mode am- 
plitude at this stage gets saturated to maintain a roughly 
constant value or register a very slow rise. The gravita- 
tional wave at the second stage falls to a very low value 



by the end of 1 year. 

Several observations in the results are in order. 

The amplitude of the r-mode first rises exponentially 
and saturates in a natural way a few hundred seconds 
after the mode instability sets in. This is unlike the case 
of [f|, where the first stage is to be terminated by hand 
and for second stage a different set of equations need 
evolving. 

For a small initial differential rotation (K ~ 0), the 
r-mode saturates at values of the order of unity. On 
the other hand, for significantly larger values of K, i.e., 
K 3> 1, the saturation amplitude can become as small as 
10~ 5 when magnetic field is absent. With the presence 
of magnetic field it can go down further to 10 -7 . Also 
the differential rotation parameter K causes to saturate 
the growth earlier than the one with K = 0. 

Higher values of magnetic field (#14=6, 10, 15 and 20, 
for example) causes the growth of mode amplitude to 
suppress drastically. The growth curves for a labelled 
with different K are merged at these higher values of 
magnetic fields (Figures Q] upto [3]) . 

The angular velocity of the star remains approximately 
constant during the first stage of evolution and abruptly 
falls to a very low value in the next few thousand seconds 
after the first stage is over. This is true for all K and 
B14. However, there exist minor differences in the an- 
gular velocity evolution curves with various K and B14 
values. Essentially, the abrupt falling of the angular ve- 
locity takes place earlier for higher K values (Figure [8] 
through HHP- Additionally, distinct trajectories are found 
at early part of the second stage in case of larger val- 
ues of K with a moderately high value of magnetic field. 
Thereafter, the curves approach each other to become 
very closely spaced, signifying a final angular velocity 
with more or less a single value at the end of the second 
stage, i.e., at the end of 1 year. 

For sufficiently high values of magnetic fields (i?i4=20 
or above), the distinct curves associated with different 
values of K all merged to a single curve. This signifies the 
fact that the effect of differential rotation on the evolution 
of r-mode instability can be erased by a large magnetic 
field. 

For investigating the gravitational wave emitted from 
the r-mode instability, we choose the same set of pa- 
rameters as in the earlier section to start the numerical 
solutions of the Equations (|36[) and (|37|) and estimate the 
gravitational wave. Figures (fT5j) through (|2"T1 depict the 
results in graphical form. 

Analysis of the graphs reveals that the amplitude of 
the gravitational wave grows exponentially within first 
few hundred seconds (first stage) to reach its maximum 
after which it falls gradually in approximately linear fash- 
ion. As the initial differential rotation parameter K 
increases, the maxima of the gravitational wave strain 
amplitude is pushed towards less time, i.e., the maxima 
reaches earlier and earlier moments as K increases. At 
the same time the magnitude of the maxima also be- 
comes lesser. As for example, the gravitational wave am- 
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FIG. 8: Angular velocity evolution for different K and S14 = 
0. 
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FIG. 9: omega vs time for different k and B14 = 10 



plitudc reaches its maximum of h max ~ 5 x 10 for 
K = and h rnax ~ lO" 30 for K = 10 12 (Figure (JTSJ) for 
the magnetic field B14 = 0. Larger values of magnetic 
field renders h max to suppress to very low magnitude. 
Also observed is the influence of magnetic field to cause 
the turning point at the maxima to be more and more 
smoothened unlike the case of no magnetic field where 
the turning point has an abrupt change. Additionally, 
larger and larger magnetic fields tend to erase the dis- 
tinction between the curves with different values of ini- 
tial differential rotation parameter K and merge all the 
curves to register a very low value of h max < 10~ for 
£>i4 = 20. The effects of magnetic fields for different 
K on gravitational wave output may be clear from the 
graphs in Figures (|19p through (f2"Tj) . Similarly the ef- 
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fects of the differential rotation factor for different mag- 
netic field strengths are depticted in Figures (fT5j) through 

It is worth mentioning here that within the model of 
evolution proposed in Ref. the frequency-domain 

gravitational wave amplitude has a spike at high frequen- 
cies, due to the fact that during the first stage of evolu- 
tion the angular velocity of the star evolves very slowly on 
the viscous timescale, leading to a quasi-monochromatic 
gravitational wave emission during the first 500 seconds 
of evolution. However, as we have seen, if one takes into 
account the influence of differential rotation, namely, the 
fact that it contributes to the physical angular momen- 
tum of the r-mode perturbation, then the angular ve- 
locity of the star evolves in the gravitational-radiation 
timescale already in the first stage of evolution. Influ- 
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ence of the magnetic field is merely to change the overall 
structure of the evolution because of the fact that the 
magnetic field term appears in the evolution equations 
([55)) and ([57)) as a separate additive term. As a conse- 
quence, the gravitational wave amplitude |/i(/)| in the 
first stage of evolution is also given by Equation ([50)) and 
we observe no spike. 
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VI. DETECTION OF GRAVITATIONAL WAVES 
FROM r-MODE INSTABILITY IN THE 
PRESENCE OF MAGNETIC FIELD 



The gravitational radiation emitted by a newly born 
neutron star, as is spins down due to the r-mode in- 
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stability, could in principle be detected by the laser in- 
tcrferometric detectors. In the following we discuss the 
detectibility of gravitational waves during the spin-down 
of neutron star by detectors like LIGO, Advance LIGO 
and VIRGO, using appropriate detecting strategies. 

In reality, neutron stars are such a complex system 
involving many other effects such as superfluidity, quark- 
gluon plasma, a very complicated magnetic phenomena, 
our knowledge of the evolution of the r-mode instabil- 
ity is insuffecient to predict the gravitational waveform 
with such an accuracy that a matched filtering signal- 
processing technique would be feasible. Therefore in this 
work, the match filtering is used only to estimate the de- 
tectibility of the gravitational- wave signal for the models 
parameters we have considered here. For this purpose, 
we construct the characteristic amplitude of the signal 
through the relation 



h c (f) = f\h(f)\ 



3.3 x 1CT 26 / 



(f/ft) (K + 2) 
13 



-1/2 



6.74 x 10 4 « 2 {f/f t y 



(52) 
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is compared with the rms strain noise in the detector, of 50< / < 1200 Hz, the curves for the noise power spec- 
tral densities of the LIGO [H, advanced LIGO ^ and 
(53) VIRGO [13] detectors are well approximated by the fol- 
lowing analytical expressions, respectively: 



,(/) = VJs^h), 



where Sh(f) is the noise power spectral density of the 
detector. For frequencies in the interval in the interval 



J 



S h (h) = Si 



where Si = 3.4 x 10" 46 Hz _1 and /i=142.0 Hz, 



h) 50 U 



where S 2 = 2.2 x 10- 47 Hz _1 , / 2 = 52.8Hz and f 3 = 421.3 Hz, and 



S h (f) = S 3 

where S 3 = 1.5 x lO^Hz" 1 , / 4 = 249.6 Hz. 



1 (h 



6 V/ 



(54) 



(55) 



(56) 
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FIG. 15: Gravitational wave amplitude vs time for different 
k and B14 = 
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FIG. 16: Gravitational wave amplitude vs time for different 
k and B 14 = 10 



In Figures (|2"2"j) through (f2"T)) the curves corresponding 
to the rms strain noises in the initial LIGO, VIRGO and 
advanced LIGO detectors are compared with the curves 
corresponding to the characteristic amplitude of the grav- 
itational wave signal for different values of differential 
rotation parameter K, the magnetic field B and for a dis- 
tance to the source of D=20Mpc. For such a distance, 
which includes the Virgo cluster of galaxies, a few super- 
novae per year is expected. 

The most striking feature of these figures is that the de- 



tcctibility of gravitational waves from the r-modc insta- 
bility of newly born neutron stars is drastically reduced 
as the initial value of differential rotation associated with 
r-modes increases. Similarly the stronger magnetic fields 
also contribute to the reduction of the signal-to-noise ra- 
tio of the these signals. (Table HI)) . 

The visual comparison in Figures (|2"2"j) through (f2"T)l 
between the charatcristic amplitude of the signal and the 
rms strain noise in the detector gives us a measure of the 
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FIG. 17: Gravitational wave amplitude vs time for different 
k and Bi 4 = 20 
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FIG. 18: Gravitational wave amplitude vs time for different 
k and Bi 4 = 100 

signal-to-noise ratio for matched filtering. A quantita- 
tive determination of the signal-to-noise ratio is obtained 
from, 

We have chosen the minimum and maximum value of 
frequency, / m in=50 Hz and / m ax=1200 Hz for estimating 
the signal-to-noise ratio here. 

We have computed the Signal-to-Noise ratios for the 
gravitational waves emitted due to the r-mode instabil- 
ity in a typical Neutron Star of mass 1.4 M Q and ra- 
dius 12.53 km. for different values of the differential 
rotation parameter and dipole magnetic field strengths 



FIG. 19: Gravitational wave amplitude vs time for different 
strengths of magnetic field and the differential rotation pa- 
rameter K — 
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FIG. 20: Gravitational wave amplitude vs time for different 
strengths of magnetic field and the differential rotation pa- 
rameter K = 10 6 

by using the standard formula (|57|) for matched filtering 
technique. The results are tabulated in the Table HH For 
the initial LIGO and VIRGO detectors, even for small 
initial differential rotation (K«0), the signal-to-noise ra- 
tio is not significant for D= 20 Mpc. However, this ratio 
can be increased if we consider sources located at smaller 
distances, but that will be at the cost of decreasing the 
number of expected supernova events per year and, hence 
the probability of a detection. 

For the advanced LIGO detector the situation im- 
proves: for small values of K, and B, the signal-to-noise 
ratio is significant even for D = 20 Mpc. Therefore, and 
since within such a distance several supernovae per year 
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FIG. 21: Gravitational wave amplitude vs time for different 
strengths of magnetic field and the differential rotation pa- 
rameter K = 10 12 
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FIG. 22: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the differential rotation pa- 
rameter K and magnetic field strength B14 = 0. 



are expected, one could hope that advanced LIGO detec- 
tors would be able to detect gravitational radiation for 
the r-mode instability of young neutron stars. 



VII. CONCLUSIONS 

We have evaluated the behaviour of r-modc instabil- 
ity including its gravitational wave detectability scenario 
when a typical neutron star goes unstable through r- 
mode, particularly in the presence of differential rotation 
and magnetic field. 

• The amplitude of the r-mode first rises exponen- 
tially and saturates in a natural way a few hundred 
seconds after the mode instability sets in. 



FIG. 23: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the differential rotation pa- 
rameter K and magnetic field strength B14 = 10. 



FIG. 24: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the differential rotation pa- 
rameter K and magnetic field strength B14 = 20. 



• Small initial differential rotation causes the r-mode 
to saturate around unity. Larger and larger differ- 
ential rotation suppresses the growth by substan- 
tial order of magnitude and also comparatively in 
shorter time scales. 

• Magnetic field plays a dominant role in dramati- 
cally suppressing the growth of r-mode amplitude. 

• In absence of either differential rotation or mag- 
netic field, the angular velocity of a typical neu- 
tron star stays constant upto a time scale of few 
hundred seconds. Subsequently it starts abruptly 
shedding in another few thousand seconds. Pres- 
ence of differential rotation causes the star to shed 
angular velocity early. On the other hand, presence 
of magnetic field causes the abrupt fall of angular 
velocity gradually smoother. 

• Initial differential rotation hardly influences the 
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FIG. 25: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the magnetic field and the 
differential rotation parameter K = 0. 



FIG. 27: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the magnetic field and the 
differential rotation parameter K = 10 12 . 
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TABLE II: Signal-to-Noise ratio of gravitional waves due to 
r-mode instability calculated for different detectors 
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FIG. 26: Characteristic gravitational wave amplitude vs fre- 
quency for different strengths of the magnetic field and the 
differential rotation parameter K = 10 . 



critical curve for r-mode instability, while the mag- 
netic field offers to shrink the instability window 
available for a typical neutron star. 

• Presence of magnetic field tends to suppress the 
effects of the differential rotation parameters. Fur- 
ther, higher strengths of magnetic field often erase 
distinctive features induced by different values of 
differential rotation parameters. 

• The gravitational wave strain amplitude emitted 
due to the r-mode instability in a typical neutron 
star grows exponentially during first few hundred 
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seconds after which it slowly falls approximately in 
a linar fashion. 

• Increased values of differential rotation cause the 
gravitational wave strain amplitude to peak at 
smaller values as well as early. Further, the mag- 
netic field smoothens this amplitude change-over 
at the peak. 



• The Signal-to-Noise ratio of the gravitational waves 
due to the models investigated in this study reveal 
that the Advanced LIGO detector with low val- 
ues of K and magnetic field upto an order of 10 15 
Gauss, will be able to detect r-mode unstable neu- 
tron stars upto a distance of 20 Megaparsecs. 
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